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Abstract: CP violating asymmetries in both the production and 
the decays of s-channel neutral Higgs bosons at colliders can 
be large. If the CP violation occurs in the muon-Higgs coupling, 
one observes a production asymmetry with transversely polarized 
beams. Likewise if the CP violation occurs in the top-Higgs 
or r-Higgs coupling, one observes an azimuthal decay asymmetry 
in n + ji— — > tt (or r + r~). CP studies at such colliders allow a 
uniquely clean way of deducing the underlying CP phases; in turn 
these parameters would significantly improve our understanding 
of baryogenesis. 
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The Large Hadron Collider (LHC) currently being planed at CERN may 
provide a reach for new physics of up to a few TeV. It is to be hoped that 
this hadronic capability will be complemented with leptonic machines. 

Due to synchrotron radiation effects it is probably impractical to con- 
struct an e + e~ storage ring (such as LEP) at these higher energies. Two 
possible leptonic colliders at the TeV energy scale that have been suggested 
are a colliding beam e + e~ machine (Next Linear Collider or NLC), which 
would operate on principles similar to the current SLC and a Muon Collider 
(MUC). Since synchrotron radiation is so much smaller at a MUC compared 
to an electron machine, a MUC has the advantage that muon beams could 
be re-circulated in a storage ring and thus muons could be reused several 
hundred times before their eventual decay. Consequently MUCs have been 
receiving increasing attention in the past few years 0-0. By the same 
token, NLC has been the subject of intense research effort B in the past 
decade and its feasibility is no longer widely questioned whereas for MUCs 
many challenges in accelerator design still need to be overcome 

For most physics applications muon collisions and electron collisions are 
virtually identical. Hence factors such as cost, and the engineering problems 
involved in building such colliders are likely to be the determining factor in 
their eventual construction. If a MUC is eventually constructed, a possible 
physics bonus may be the ability to study neutral Higgs bosons in the s- 
channel 0, |j . 

In this paper, we shall assume that at least one neutral Higgs scalar, which 
we will denote 7i, will be discovered in the mass range 100 < < lOOOGeV^ 
and that the muon collider will be tuned to study such a resonance in the 
s-channel so that y/s = m^. In this way the feeble coupling of the /i to the 
Higgs may be partially compensated by the resonance enhancement and so 
the signal of an s-channel Higgs boson may be studied in detail. Bearing this 
in mind we consider how the CP phase of its couplings to fermions may be 
investigated. We find that the method gives a uniquely clean way of deducing 
CP violating phases in the coupling of the neutral Higgs to leptons and to 
the top quark. 

As is well known, at the present time, the only form of CP violation 
experimentally observed is in the K° system. It is thought that these results 
can be explained through the standard model's phase in the CKM matrix [^], 
a hypothesis which will be tested at B factories now being built. However, 
strong theoretical arguments exist that suggest that additional forms of CP 
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violation are needed to explain the baryon asymmetry of the universe [H]. 
Indeed one class of models which has been shown to give reasonable baryon 
asymmetries are models involving CP violation in the Higgs sector. 

Motivated by these considerations we shall investigate two methods of 
observing a CP violating signal in such a Higgs resonance at a MUC. First, 
if the /i beams are transversely polarized at different inclinations, there is 
a CP violating production asymmetry proportional to sin0 M , <f)^ being the 
azimuthal angle between the polarization directions. This is sensitive to CP 
violation in the Higgs-/i/i coupling. The second observable we consider is 
an analogous asymmetry in fi + fi~ — > tt (or t + t~). In the latter case, for 
instance, the signal is proportional to sin0 t where <pt is the azimuthal angle 
between the polarization of the top quarks; <pt may easily be determined 
if the tops decay leptonically. Polarization of the colliding muons is not a 
necessity for this study, although longitudinal polarization of the beam can 
be used to reduce the background. In this case one is sensitive to CP violation 
in the Higgs-tt (or t + t~) coupling which relates directly to some proposed 
mechanisms of baryogenesis in the early universe || . 

In all three cases mentioned in the preceding paragraph the asymmetries 
may be large. However, to a certain extent, they will be obscured by standard 
model backgrounds, as we will discuss below. We find that observation of 
these asymmetries will certainly be possible if the Higgs is relatively light or 
narrow leading to a relatively large production rate. Furthermore, in many 
theoretical scenarios, even if the Higgs is not that light, it can have a large 
production rate and it can lead to appreciable asymmetries. 

Although the underlying model we have in mind may contain several 
Higgs doublets we would like to cast our analysis in a relatively model inde- 
pendent way. Let us assume that a single Higgs is under experimental study 
which we will denote 7i. We will parameterize its coupling to a fermion / by 

Cms = c° ff xs^ (i) 

where C ^ is the coupling in the standard model (with one Higgs doublet). 
In such theories with an extended Higgs sector the field TC may either be 
a scalar if or a pseudoscalar A. If TC = A then it does not couple to two 
gauge bosons while if TC = H we will parameterize the coupling to two vector 
bosons (VV = ZZ or WW) as: 

C H vv = Cy V cos a. (2) 
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where Cy V is the coupling of the SM Higgs to the gauge bosons and a is the 
angle between the observed Higgs TC and the orientation of the vacuum in 
the Higgs space. 

To serve as an illustration we will consider three sample models for the 
TC in question: 

• 1: TC — H with a = A/ = tt/4 and Xf = 1 f° r & U fermions. 

• 2: TL = A with A/ = it /4 and Xf — 1 f° r all fermions. 

• 3: 7i = A with A/ = 7r/4 and x« = Xd = 5 and Xu = 1/5. 

The last case is motivated by models of the type that, for instance, are 
used for the Higgs sector for supersymmetric (SUSY) scenarios ||. In such 
a case one might want to partially explain the large top mass by making 
the corresponding vacuum expectation value (vev) large giving rise to xi = 
Xd — v%/vi — 1/Xu for large Vijv\. Although minimal SUSY models do not 
have CP violation in the Higgs sector, more complicated Higgs sectors or 



non- minimal SUSY models do readily admit CP violation [10 . 

It is difficult to observe the effects of the Higgs resonance except if the 
accelerator is precisely on shell 0, H. If s is the center of mass energy and 
s = then defining = Br(7Y — > and an = cr(yU + /i _ — > 7i), the 

cross section for Higgs production is: 

An 

a n = —2 Bp. (3) 
It is useful to compare this with er = <j(/i + / u _ — > 7* — ► e + e~) so we define 

R{H) = ^ = ^ 2 B, (4) 

CT "e 

where a e is the electromagnetic coupling. It is clear in the above that the 
key to having a large rate for Higgs production is for its width to be as small 
as possible (thereby increasing B^). The main possible decay modes of the 
Higgs which we consider are H — > ZZ, H — > W + W~, H — > ti, TC — » bb, and 
TC — > r + r _ . The last two become important for the case when < 2mw 
The decay rates to these modes, given the above couplings, can be readily 
calculated by using the results that exist in the literature P|: 
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r(W-tt) = 39 £ mt T H (3M + (l-ffism\ t \ X 2 t K t 



F(H ^ bb) = -^^XIK* 



nn^WW) = + (5, 

where /3j = y 1 — Amf/m-H- The and VTW 7 modes are only relevant if 
7i = H . In these equations i^, i^t are QCD corrections given, for example, 
in Ref. @. 

In comparison, the width of the Higgs decay to is 

In Figure 1 we present the value of R(7i) at s = assuming that a = 45° 
for the H = H and H = A function of m w . 

In deriving equation (|) above we assumed that the energy of the collider 
could be controlled to a precision much finer than the width of the Higgs. 
While this should be generally true if m H ^> 2m w or 2m t , below this thresh- 
old it may well not be the case. To take this into account let us assume that 
the actual value of s is uniformly distributed in the range: 



m 



H 



1-5) < s <m 2 n (l+5) (7) 



The beam spread thus defined will lead to an observed rate of Higgs produc- 
tion described by: 



R{H) 



■ arctan 



H 



R{H) (8) 



In Figure 1 we also show R{7i) using the value of 5 = 10~ 3 . 

Let us now consider the production of — ■> H with transversely 
polarized \x beams. We take the z-axis in the center of mass frame to be the 
direction of the pT beam to be polarized in the x direction. Assume further 
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that the polarization of the /i + beam is inclined at an angle of M . The cross 
section in this case becomes: 

cr(0 M ) = (1 - cos 2A M cos M + sin 2A M sin M )er o (9) 

where do is the unpolarized cross section. Then the CP odd asymmetry in 
the production will be given by: 

ct(90°) -<t(-90°) x , , 

A u = —. ( f = sin2A H 10 

a(90°) + a(-90°) V ; 

Thus the underlying CP violating phase A^ becomes cleanly measurable. 
Since this phase (A At )is unconstrained, large asymmetries approaching 100% 
are possible. As we shall discuss below, this asymmetry will be diluted 
somewhat by the standard model backgrounds. 

Let us now turn our attention to the analogous type of asymmetry in the 
decays 7i — > ti or t + t~ ||11|| . In order to observe a CP violating correlation 
in the spins of the top quarks or r-leptons, we clearly need to look at the 
correlations among their decay products. 

Consider first the determination of the spin of a fermion in general. Sup- 
pose / decays to XY where X, Y may be a single particle or a multiparticle 
state. Let e x be the "analyzing power", i.e., the degree to which the mo- 
mentum of X is correlated with the spin of /, defined by 



-x 



3 < cos9 x > ■ (11) 



where 8x is the angle between the spin of / and Px, in the / rest frame. 

We discuss the decay of the top quark first. If the top decays semi- 
leptonically, t — > bW — > bl + vi (for I = e, //) then ef = 1 [I2~|. For hadronic 



decays of the top one could obtain e = 1 if one could tag which jet was the 
d type (thus playing the role of the positron). Experimentally this is hard 
to do. We will therefore only use the inclusive W momentum to determine 
the polarization for the hadronic decays. The analyzing power in this case is 
thus: 



' W m\ + 2ml 



^0.39 (12) 

Next consider the case of the r lepton. If the r decays semi-leptonically, 
then ej — — | |13| (note that if one could observe the V\ then e = 1). In JL3| 
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it is shown that the inclusive hadronic decays of the r can give = — . 42 
(h stands for inclusive hadron states). If one analyses the detailed structure 
of the hadronic decays one can obtain better results [Hj] but we will just 
use this number in our discussions to follow. 

Let us define a coordinate system in the Higgs center of mass frame where 
the z' axis is in the direction of the / (i.e., t or r) momentum. Let us now 
consider the / decays via / — > X^Yi and the / decays / — > XjYj. We define 
the angle 0y to be the azimuthal angle between the pxi and the p X j projected 
into the x' — y' plane: 

. ( , v Pxi x p Xj ■ p f 

sin ^ = p 13 

\PXi\ \Pxj\ \Pf\ 

Taking into account e{ , the differential distribution in (fiy is (note that 
ef = -ef): 



dT IT 2 t t 7T 2 f f 

r^JT~ = 1 + j^titjPf cos2\ f cos (pij + — efe^/ sin2A/sin0„ (14) 



J '.) 



In eqn. (|Tj), for f — t (so that the threshold mass effects must be retained), 
p and i] are given by: 



Pt 



l-& 2 -(l + /f)cos2A t 
cos2A*(l + A 2 - (l-/3 t 2 )cos2A t ) 



Vt ~ l-(l-A 2 )cos 2 A, (15) 

On the other hand, if / = r then the above simplifies to p T — rj T — 1. The 
resultant CP violating asymmetry which we will consider is: 

f = r(sin0j 3 - > 0) - r(sin0jj < 0) 

ij r(sin0 ij > 0) + r(sin0 ij < 0) [ ) 
Thus using the above distribution we find: 

A ii = l4^Vf sin2X f ( 17 ) 



Of course we would like to optimally combine the resultant asymmetries 
from different pairs of modes {ij}- In order to do this we would like to choose 
weights Wij normalized by: 
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», ; Y^is.Bj (is) 

where Bi is the branching ratio of / — ► i and the summation is over all the 
modes under observation. Thus the total asymmetry 

A f Y^if.jM^lhBj (19) 

is to be maximized. This is done by taking Wij oc A{j [16|. Applying this to 
the asymmetries defined in fll6|) we find that the total asymmetry is given by 

Af = ^(Sf)%sin2\ f (20) 

where 



^ = VE^) a ( 21 ) 

Selecting the top quark and the r lepton modes mentioned above we obtain 

Jt = 2 
M 9 



numerically (using m t = 17QGeV : BJ = .36, = .64, B l e + 5* = | and 



5* +5* 



5* = .58 £ T = .39 (22) 

In the case of top decay, if the quark or the anti-quark decay semi- 
leptonically, the 4-momenta of the undetected neutrino (s) may be recon- 
structed by using the mass constraints of the top and the W. For the decay 
of the r lepton, it is less obvious that this can be done. For instance, if 
one or both of the r leptons decays leptonically, there are at least 3 z/'s in 
the final state and the kinematics does not provide sufficient restrictions to 
reconstruct the event. Even if both r's decay hadronically, although it would 
appear that there is enough information to reconstruct the event, in fact 
solving the kinematics gives rise to a "quadratic ambiguity" resulting in an 
ambiguity in the sign of sin 0y . Therefore with this information alone the 
asymmetry A? cannot be evaluated. 

These problems with the r may be resolved as follows given a reasonable 
vertex detection capability. 

Suppose Xi is a particle emerging from the r~ and Xj is a particle emerg- 
ing from the r + . In general the tracks of these two particles will not intersect. 
Let us define the vector V{j to be the displacement from the track of Xj to the 
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track of X{ at either their closest approach or the closest approach inferred 
from extrapolating the portion of the tracks which may be observed. Typi- 
cally the magnitude of Vij will be of the order of 2cr T = 180/xm (in spite of the 
fact that the displacement of the decaying r's is of order 2cr T 7 = O(lcm)). 
We then define the angle ^ij to be the difference in the azimuthal angle 
between p x % and p X j with respect to an axis in the direction. 

Recall that the asymmetry Ay is just the expectation value of cr(sin(0 ij -)) 
where: 



a ( x ) = \ 1 if „ ^ n ( 23 ) 




A vertex detector may allow the determination of <r(sin0) in a r + r event 
using one or more of the following four methods: 

• 1: Since the displacement of the r decay from the interaction point is 
on the order of cr T 7 = 0{lcm) it is not unreasonable to assume that, 
depending on detector design, in many events either the r + or the r~ 
will form a track in the vertex detector. If this is the case then the z' 
axis is determined and <pij is obtained directly from the momenta of Xj 
and Xj. 

• 2: If Xj or Xj consist of multiple charged tracks, running these tracks 
back to their common decay vertex again determines the z' axis if the 
machine is configured so that the interaction region is in a single point. 
More generally if the interaction region is not a well defined point, one 
may use the decay vertex and the combined momentum of the multiple 
tracks to define the trajectory of the center of gravity of the system 
and use method (3) below treating this trajectory as an effective single 
track. 

• 3: If both Xi and Xj are single charged tracks (e.g., r~ — > pTv^ or 

— > 7r~z/ T ) then 



a (sin faj) = a (sin faj) (24) 

so that the exact decay vertex need not be identified. To see why (|4j) 
should be true, let x\ be any point on the track of Xi and X2 be any 
point on the track of Xj. If u = x\ — x~2 then u-(piX pj) = v%j ■ (pi xpj)- 
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If the t~ decays at time t\ and the r + at time t 2 then the vertex points 
Vi, Vi are V\ = m~Hip T - and V 2 = m~ 1 t 2 p r + therefore V% — V 2 = 
m^iti + t 2 )p T - for the r~ and r + respectively. Since the point V\ lies 
on the Xi track and the point V 2 lies on the Xj track, 



a sin 1 



• Pi x p,) 

(7((Vl - V 2 ) - Pi X p}) 

o-(p7 • pi x p,) 

er(sin (25) 



• 4: If both r's decay hadronically we can use the method of [17]. In 
this instance the hadronic states will contain charged tracks Xj — > 
! . . . X i>n and Xj — > X^x . . . X J)Tn . For any pair of charged tracks 
{Xj fc, X^y} we define 0^ as above. In this case we first kinematically 
solve for the neutrino momenta up to the "quadratic ambiguity" men- 
tioned above. Note that this is in fact a parity ambiguity so that only 
one of the reconstructions will give the correct values for sin 0j ^ ^y. The 
false reconstruction will give exactly the wrong sign for each of these 
angles providing multiple checks on the parity of the reconstruction. 

The only decay combination of modes which is not obviously susceptible 
to methods 2-4 (although method 1 will work in all cases given a vertex 
detector sufficiently close to the beam axis) is Xj = 7r~7r° and Xj = l + 
(or vice versa). Let us assume that the 2tt mode is dominated by the p 
resonance. Since, using equation ( |T2"D generalized to r — >• pv gives e T p = .45 
so that e T p e[ ~ .15 ~ (£ T ) 2 . Thus totally ignoring these modes is roughly 
equivalent to reducing the cross section by 2B T p BJ w 10%. 

The above asymmetries must be observed against the CP even standard 
model backgrounds p + p~ — > bb, ti, ZZ and W + W~ . Therefore, using the 



expressions for the rates for these reactions given in Ref . [18j , we numerically 
study their impact on the CP violating signals that we are seeking. 

When looking for decay asymmetries of the Higgs the ability to control 
the longitudinal polarization of the beams may help to increase the signal 
over the background |19[. This is because the Higgs production flips the 
chirality of the muon while all standard model interactions preserve it. Thus, 
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for instance, if both beams are left polarized with polarization P then the 
Higgs production is enhanced by 1 + P 2 while the background is reduced by 
1 — P 2 [[TJ|. In the numerical results below we will illustrate the effects of 
polarization for a few cases. 

Let R® be the standard model contribution to a specific final state j. Let 
Lq be the integrated luminosity per year of a given accelerator. For a final 
state / let us denote by yf^ the number of years needed to accumulate a 
3-0" signal for the CP asymmetry. This will be given by 

y M = 9 R° + RnBr(H^j) 

We will use the notation yf a for the CP asymmetry in the production ob- 
tained by monitoring the decay of the Higgs to the final state % and y 3a for 
the azimuthal asymmetries described above by observing the decay of the 
state j (i.e., j = tt or tt). 

In the following we will take the luminosity to be 10 34 cm _2 s~ 1 so as- 
suming that a year is 10 7 s (taking 1/3 efficiency), Co = 10 41 cm~ 2 yr~ 1 . 
Numerical results for the three cases of Higgs couplings given below eqn. (Q) 
are summarized in Figures 1-4. For definiteness in all the cases we will take 
X f = it/A. 

In Figure 1 we consider the overall Higgs production cross sections in 
relation to fi + fi~ — > 7* — > e + e~ . In each of the three cases we show Ry, and 
Rti with 5 = 10 -3 . Note the pronounced effect of the vev ratio in case (3) 
due to the fact that the decay to tt is suppressed and the bb and modes 
are enhanced. 

In Figure 2 we consider y 3cr for production asymmetries as monitored 
by various final states in case (1). Below the WW threshold y 3cr is in the 
10~ 3 — 10~ 4 range so observation should be relatively easy. Indeed, in this 
instance, even a machine luminosity of 10 31 -10 32 cm _2 s _1 may be adequate. 
Above this threshold y\y W rapidly increases and passes through 1 at m n 
250 GeV. 

In Figure 3 we consider the same quantities for case (2). Here we take 
5 = 10" 3 in all the curves. In this case the WW and ZZ thresholds are 
absent so a small value of y 3cr persists up to the tt threshold. So again for 
this mass range the production asymmetry is a very sensitive probe and 
observation may be feasible even with a machine of £ ~ 10 32 cm~ 2 s -1 . After 
the it threshold y 3 f rapidly grows passing through 1 around ~ 400 
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and observation starts to become rather difficult. We also show the decay 
asymmetries y 3 £ and y^f. y 3 £ is about 10 below the ti threshold, yff grows 
rapidly from 1 to 100 above threshold. The effect of polarization is also 
illustrated in this figure for these two instances of y^ T and y^f. The lower 
dotted and dashed dotted curves are with P = 0.9 (the upper ones with 
P = 0). We thus see that polarization can improve the effectiveness of the 
decay asymmetries by about one order of magnitude. 

In Figure 4 we plot the same quantities for case (3). In this case y 3a 
ranges between 10~ 5 and 10~ 3 over the entire mass range so the production 
asymmetries are excellent probes for < 800 GeV. Indeed for this case 
the production asymmetries appear detectable even with C ~ 10 31 cm _2 s _1 . 
The decay asymmetry yff is around 10~ 2 and y 3 " is in the range .1-1 below 
the ti threshold and 1-100 above it. We have also shown the effect of using 
polarization P = .9 on yf^ and ytf a \ Again polarization seems to improve 
the situation by about an order of magnitude. Thus the decay asymmetry 
in the rr channel is quite effective for < 300 GeV; for a heavier Ti the ti 
channel seems very promising and C ~ 10 32 cm _2 s _1 may well be adequate. 

Summarizing, this study shows that production asymmetries represent a 
promising way to study CP violation in the couplings of Ti in all the three 
cases for at least some range of Higgs masses. Decay asymmetries appear 
less effective and probably are most useful in models where there is a large 
ratio between the vacuum expectation values. Longitudinal polarization of 
the muon beams helps to significantly reduce the backgrounds to the decay 
asymmetries and can render them quite viable in a variety of Higgs models. 

This work was supported by US Department of Energy contracts DE- 
AC03-765F00515 (SLAC) and DE-AC02-76CH0016 (BNL). 
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Figure Captions 

Figure 1: 

The value for R-^ as a function of for the following three cases: (1) 
U = H, a = tt/4 and Xf = 1 (solid) (2) H = A and x/ = 1 (dashes) (3) 
H = A and xi — Xd = 5 and x„ = 1/5 (dots). In each case the upper branch 
represents the result for y/s = while the lower branch is the result with 
an energy spread given by 5 = 10~ 3 . 

Figure 2: 

The value of for case 1 (see caption to Figure 1) assuming a = tt/4 = 
A M . y^ is shown with the solid line, y^w with the dot-dashed line, y^z 
with the dotted line, and y tl ° with the dashed line. (Note y 3a is the number 
of years needed to accumulate a 3a production asymmetry. Note also that 
the horizontal line, y 3a — 1, is drawn to serve as a point of reference.) 

Figure 3: 

The value of y^ for case 2 assuming A M = A T = X t = 7r/4. y^ is shown 
with the solid line, and y[^ with the dashed line. The results for decay 
asymmetries are also shown: y^ with the dot-dashed line and y^ with 
the dotted line. In both of the latter cases the lower curve is for P = .9 and 
the upper curve is for P = 0. See captions to Figures 1 and 2. 

Figure 4: 

The value of y^ for case 3 assuming A M = A T = At = tc/4. y^ is 
shown with the solid line, and y^ with the dashed line. The results for 
decay asymmetries are also shown: y^ with the dot-dashed line and y^ 
with the dotted line. Again in both of the latter cases the lower curve is for 
P = .9 and the upper curve is for P = 0. See captions to Figures 1 and 2. 
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